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We study one-dimensional p-wave superconductors capacitively coupled to a microwave stripline 
cavity. By probing the light exiting from the cavity, one can reveal the electronic susceptibility of the 
p- wave superconductor. We analyze two superconducting systems: the prototypical Kitaev chain, 
and a topological semiconducting wire. For both systems, we show that the photonic measurements, 
via the electronic susceptibility, allows us to determine the topological phase transition point, the 
emergence of the Majorana fermions, and the parity of their ground state. We show that all these 
effects, which are absent in effective theories that take into account the coupling of light to Majorana 
fermions only, are due to the interplay between the Majorana fermions and the bulk states of the 
superconductors. 
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I. INTRODUCTION 


Condensed matter systems are an endless resource of 
emergent physical phenomena and associated quasipar¬ 
ticles. Majorana fermions, which are particles that are 
their own antiparticles and which have been first pro¬ 
posed as particles in the context of high energy physics, 
emerge beautifully as zero energy excitations in con¬ 
densed matter setups mm- Specifically, they are pre¬ 
dicted to occur as zero energy excitations in solid-state 
systems, such as genuine p-wave superconductors [3H5], 
or engineered from topological insulators [6j, semiconduc¬ 
tor wires in a magnetic field EH], or in chains of mag¬ 
netic atoms [IMS], all in the proximity of s-wave super¬ 
conductors. These exotic objects are robust against lo¬ 
cal perturbations and, moreover, they obey non-Abelian 
statistics mmm under braiding operations, thus rec¬ 
ommending them as qubits for the implementation of 
topological quantum computation. 

Electronic transport is the foremost experimental tool 
for investigating the Majorana fermions physics but al¬ 
ternative, non-invasive , methods that preserve the quan¬ 
tum states would be highly desired to address these ob¬ 
jects. Cavity quantum electrodynamics (cavity QED) 
has been established as an extremely versatile tool to ad¬ 
dress equilibrium and out-of-equilibrium electronic and 
spin systems non-invasively [EM]. Majorana fermions, 
too, have been recently under theoretical scrutiny in the 
context of cavity QED physics [281437] . However, most 
of these studies dealt with effective low energy models 
that involved Majorana fermions only, leaving the bulk 
physics, which is at the heart of the Majorana physics, 
largely unexplored. 

The basic idea behind cavity QED with electronic sys¬ 
tem is that it allows one to extract various properties of 
the latter, such as its spectrum and its electronic distri¬ 
bution function, from photonic transport measurements, 
as opposed to electronic transport. Such photonic trans¬ 
port is quantified by the complex transmission coefficient 
r = Aexp(i</>) that relates the output and input photonic 
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FIG. 1. A sketch of the system: a one dimensional system (red 
rectangle) is placed at the maximum of the electrical field 
(green straight arrows) inside a superconducting microwave 
cavity (blue). The electromagnetic field inside the cavity is 
probed by sending input fields of amplitude and phase Ai n 
and <4>in , respectively, and measuring the field at the end with 
A out and front- The difference between the two gives a direct 
access to the electronic correlation function in the wire (see 
text). The presence of Majorana end modes in the finite wire 
(black curves) is also signaled in the cavity response. 


fields as depicted in Fig. [l] In the weakly coupled limit, 
one finds [58} [39 ] . Appendix |B| 

t(cj) = - 7 - r —- , ( 1 ) 

— i(uj — uo c ) + k — i n(cj) 

where uj c and k are the frequency and the escape rate 
of the cavity, respectively, while n( cj) is an electronic 
correlation function that depends on the actual coupling 
between the two systems, and which contains informa¬ 
tion about the spectrum of the electronic system. The 
phase and amplitude response of the cavity close to res¬ 
onance uj ~ uj c are related to the susceptibility n(cj) as 
follows: 5fr = U'(uj)/k and 5A/Ai n = TV"(w)/n, where 
5(f) — (front (frim 5 A — A.i n A OU f , and n (cj) — 7Ze[U(uj)\ 
(n"( uj) = Imp(w)]) is the real (imaginary) part of the 
susceptibility. 

In this paper, we evaluate the function H(uj) first for 
the simple case of a one-dimensional (ID) p -wave super¬ 
conductor described by the Kitaev chain and then for 
more realistic model of a ID topological semiconducting 
wire in proximity of a superconductor. We assume in 
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both cases that these ID systems are coupled to a mi¬ 
crowave cavity, as showed schematically in Fig. [l] We 
address various physical situations for this coupling and 
show that such a method allows us to ascertain the topo¬ 
logical phase transition point, the occurrence of Majo- 
rana fermions, and the parity of the ground state, all 
in a global and non-invasive fashion. The paper is struc¬ 
tured as follows. In Sec. [TTJ we describe our model Hamil¬ 
tonians for the two systems under consideration and dis¬ 
cuss the coupling between the microwave photons and the 
electrons in the ID topological systems. In Sec. m we 
show how the optical transmission through the cavity is 
able to probe the topological phase transition. In|TV| we 
demonstrate that the cavity allows to detect the occur¬ 
rence of Major ana fermions and the parity of the Majo- 
rana fermionic state in a non-invasive fashion. Finally, in 
Sec.[V]we provide a brief summary of our results. Techni¬ 
cal details of the calculations are given in the appendices. 

II. MODEL HAMILTONIAN 

In the following we will consider various models of p- 
wave superconductors coupled to a microwave (supercon¬ 
ducting) cavity, such as the Kitaev p-wave superconduc¬ 
tor model, and the spin-orbit coupled nanowire subjected 
to a magnetic field and in the proximity of an 8-wave su¬ 
perconductor. 

The general Hamiltonian for the one-dimensional sys¬ 
tems we consider here is of the form: 

H-sys — H e i T H e i_ c + Hph , (2) 

being the sum of the electronic Hamiltonian, its capac¬ 
itive coupling to the cavity, and the free photon field, 
respectively. While the electronic term is model specific, 
and it will be discussed below, the last two terms read: 

N 

U e i- C = a ^2 hi (a + a f ), (3) 

i=1 

and 

H ph = w c a) a . (4) 

In Eq. §, ^ (a) is the photon creation (annihilation) 
operator, respectively, a is the electron-photon coupling 
constant that couples to the charge density h. This 
merely acts as to shift the chemical potential. In Eq. 0, 
uj c is the frequency of the photonic mode (setting h = 1 
throughout). Such a model could be realized experimen¬ 
tally by coupling a spin-orbit nanowire in the presence of 
a Zeeman field to a nearby 8-wave superconductor mini- 
In the present setup, which is based on a microwave su¬ 
perconducting stripline cavity, the 8-wave superconduc¬ 
tor that induces superconducting correlations in the wire 
could be a part of the underlaying cavity. For exam¬ 
ple, the nanowire could be tunnel-coupled to the central 
superconducting material showed in Fig. [l] We have a 


considered a global capacitive coupling between the elec¬ 
tronic nanosystem and the cavity electric field. Such a 
coupling can be justified by a full microscopic approach 
(see Appendix |A| for details and also Ref. m that pro¬ 
vides a microscopic description of the electric coupling 
between electrons in a nanocircuit and cavity photons). 

By solving the equation of motion da/dt = —i[a, H sys \ 
for the photonic field iteratively up to second order in 
a with respect to the cavity frequency uj c [39], we find 
for the correlation function n(cj) in Eq. 0 in the time 
domain 

n (t — t') = — ia 2 0(t — t')([hi(t),hi(t ')]), (5) 

being the total charge susceptibility of the electronic sys¬ 
tem (which can be here a ID p —wave superconductor 
or a topological ID wire). In Eq. ([5|, we introduced 
hi(t) = W(t)hU(t), with n being the total number of 
electrons operator and U(t) = exp (—iH e it) the evolu¬ 
tion operator for the electronic system. We assume zero 
temperature limit (T = 0) so that the average (...) is 
taken over the superconducting ground state. Note that 
n(o;) = dt exp (ioot)U(t) and that T1(uj) = 0 in the 
absence of superconductivity (A = 0), i.e. there are no 
effects from such a coupling for a wire in the normal 
state. We detail below the models for both topological 
ID systems we consider in this paper. 

A. Kitaev chain 

The simplest model of a p-wave superconductor that 
hosts Majorana fermions is the Kitaev chain |T. There¬ 
fore, we first consider for the electronic part in Eq. (J2|, 
the Kitaev Hamiltonian that reads: 

N i N—l 

H ei =-T c i c i“oE (t c ] c i+1 + A c j c j+l + h.c.), 

3 = 1 1 3=1 

(6) 

where t is the hopping parameter, A is the p-wave super¬ 
conducting pairing potential, p is the chemical potential, 
and N is the total number of sites. Also, cj(cj) is the 
creation (annihilation) electronic operator at the site j. 
Note that the electronic operators are spinless, and the 
electronic density is given by h = XlyLi c ] c j- In the 
present setup, which is based on a microwave supercon¬ 
ducting stripline cavity, the 8-wave superconductor that 
induces superconducting correlations in the wire could 
be a part of the underlaying cavity. For example, the 
nanowire could be tunnel-coupled to the central (super¬ 
conductor showed in Fig. [l] The fact that microwave 
photons effectively couple only to electrons of the Kitaev 
chain is accounted for in Appendix |A| 

The Kitaev Hamiltonian in Eq. ([6^ can be easily diag¬ 
onalized. The susceptibility in this case can be found 
by simply substituting the expression for the density 
in Eq. 0 with the one corresponding to the Kitaev 
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model. We will discuss its physical content in Sec. Ill 
and Sec. EO 


B. Spin-orbit coupled nanowire 


A realistic system that can emulate, in some limits, 
the Kitaev chain consists of a nanowire with a spin-orbit 
interaction, subjected to an external magnetic field, and 
coupled by proximity effect to an 5-wave superconduc¬ 
tor [ZHU The entire system is then assumed to be (ca- 
pacitively) coupled to the microwave cavity. The tight- 
binding Hamiltonian for the nanowire with spin- 
orbit (SO) interaction in the presence of the magnetic 
field reads m 


N- 1 N 

H e( = “* E C i+1 JclP c 3P ~ 

3 = 1 3=1 

N N—l 

+ A E c )Ai - *7 E c j+i 

3=1 3=1 

N 

-^Eto+ Lc - (?) 

3 = 1 

where, as before, t and g are the hopping amplitude and 
the chemical potential, respectively, 7 is the spin-flip hop¬ 
ping amplitude (or the spin-orbit coupling), A is the s- 
wave pairing potential induced by proximity, Vz is the 
Zeeman splitting energy (Vz = —ggsB/ 2 , with g and B 
being the ^-factor and external magnetic fields, respec¬ 
tively). Also, Cja- (cj a ) are the annihilation (creation) 
operators for electrons at site j and spin a and 

cq, with i = x,y,z are the Pauli matrices that act in 
the spin space. This model accounts thus for spinfull 
electrons. Note that we assumed the spin-orbit field and 
the magnetic field to be orthogonal. The coupling to the 
cavity is again capacitive, and the density reads in this 
case hi = J2 a c io- c ^- In order to find the susceptibility, 
we need to substitute this expression for the electronic 
density in Eq. ([5|, and we will discuss the various cases 
in the following sections. 


III. TOPOLOGICAL PHASE TRANSITION 


Next we will show that the topological phase transition 
can be inferred from the cavity response from the trans¬ 
mission t(cj) via the susceptibility n(cj). This function 
can be calculated straightforwardly in the case of a closed 
ring, i.e. for periodic boundary conditions (PBCs), so 
that c/v+i = ci for the Kitaev chain (c/v+ i a = c\ a for 
the SO nanowire). 


A. Kitaev chain 


For PBCs, we can utilize the Fourier description for 
the electronic operators: Cj = 1/V^Efe e * fej ' c fe, with 
k = 2i:n/N (assuming the lattice spacing d = 1 thereon), 
with n = 1... N. For more details see Appendix O By 
doing so, we can readily write down the electronic Hamil- 
tonian H% = H% dG (k), with 

H Bdo( k ) = (“t cos k - n) 7 - A sin fc 7 , (8) 

where r k = (r k , r k ,r k ) are Pauli matrices that act in 
the Nambu (particle-hole) space, i.e. on the vectors &k = 
(c/c, cl fe ). The coupling to the cavity, on the other hand, 
simply reads 


Hei-c — g T T Z ( flt + a ) > (9) 

k 

so that the susceptibility in the time domain can be writ¬ 
ten as: 

n (t) = —ia 2 y^(0|[r*(i),7(0)]|0), (10) 

k 


with r k (t) = e lH BdG( k ) t T k e -i H BdG( k ) t , Utilizing this de¬ 
scription, after some lengthy but straightforward calcu¬ 
lations, we obtain for the susceptibility (in the uo space): 


n(w) = -a 2 ]T 

/c>0;p=± 


(A sin k) 


p 


El 


id + 2 pEk + ig 


( 11 ) 


where E% = yj (—t cos k — g) 2 + (A sin k ) 2 is the Bogoli- 
ubov spectrum of the ID p-wave superconductor [from 
diagonalizing the BdG Hamiltonian in Eq. |8|] [T: and 
g is a small positive number that accounts for causality. 
For t = A, the imaginary part U"(uj) acquires a simple 
analytical form, and it is given by 


n» = 


a 2 tN , 

2fiud 


\ 


1 - E 


(ud/2) 2 - t 2 - /i 2 


2 .,2 


4t> 


( 12 ) 


for |t + ji\ < ud/2 < |t — /jl\ and is zero otherwise. The 
topological phase transition takes place at \g\ = t, with 
the system being in the topological (trivial) phase for 
\/jl\ < t (|/i| > t). In Fig. [2] we plot U"(uj) (main plot) 
and n'(u;) (inset) as a function of the chemical potential g 
for various values of the cavity frequency uo. We see that 
this function shows a large peak at the transition point 
(|/i| = t), which becomes narrower and more pronounced 
for smaller uj (compared to the gap A). Physically, this is 
due to the fact that the electronic levels close to the zero 
energy have larger curvatures, i.e. they are more suscep¬ 
tible close to the phase transition point. The real part 
also serves for detecting the phase transition, although 
not as directly as the imaginary part, as shown in Fig. [2j 
where the phase transitions are inferred from the kinks 
in this function. We have checked that the same peak 
structure holds for the cases when A /t, too, the only 
modification being a shift in the scale for cj, which should 
be of the order of cj ^ A. 
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FIG. 2. The imaginary part of the density-density correlation 
function [hGfh;)] as a function of /x for the Kitaev model. The 
topological phase transition takes place at /x — — 1, where this 
function reaches its maximum, indicating the transition point. 
Inset: The real part of the density-density correlation function 
[II^cj)], which also shows features (kinks) around the topo¬ 
logical phase transition point. The full, dashed, dot-dashed 
curves correspond to the uj — 0.2, 0.3, and 0.4, respectively. 
We assumed t = A = l,iV = 50 and all energies are expressed 
in terms of t. 



Vz[meV] 


FIG. 3. The imaginary part of the density-density correlation 
function [II // (a;)] as a function of the Zeeman splitting Vz 
for a SO coupled nanowire in case of PBCs for parameters 
L — 2 /im, t — 0.5 • 10 -2 eV, a — 0.4 meV, p — —10 -2 
eV, A — 0.25 meV, N — 80, and uj =s 0.1 meV (cj = 0.2 
meV in the Inset). The topological transition takes place 
for Vz ~ 0.25 meV, for which the susceptibility reaches its 
maximum. The emergence of a second peak is due to the 
resonance condition around the external gaps (k ~ /cf). 


B. Spin-orbit (SO) coupled nanowire 


The case of a realistic SO coupled nanowire is more 
complicated that the Kitaev model showed above, and 
so is the evaluation of susceptibility. This is so be¬ 
cause the SO coupled wire has four bands (because of the 
spin), instead of two, and a more complicated quasipar¬ 
ticle spectrum. Nevertheless, writing the electronic op¬ 
erators in the Fourier space as c Jcr = 1/y/N e %/ej c/ co -, 

we can write again the electronic Hamiltonian as = 
with 


H^g = K -t cos k - n) + 7 sin k <t z ]t* + Vz<J x + At* , 

(13) 

and the coupling to the cavity the same as in Eq. 
However, the expression for U(t) becomes rather cum¬ 
bersome for the general case and to get some analytical 
insights we need to resort to approximations. For that, 
the Hamiltonian can be put in a different form by the use 
of a unitary transformation (see Appendix |D|): 


H w 


BdG(ty — 
+ 


—t cos k — /i + y (7 sin k) 2 + Vj cr z 
Ay sin k k AVz 


^/(qsin k) 2 + Vj x ^/(qsin k) 2 + V| v y ’ 

(14) 


while the H e i_ c stays unchanged. Progress can be made 
if we assume the limit of large magnetic field, Vz A, /x, 
in which case we can neglect the last term in the above 
Hamiltonian. By doing so, we recover two copies of the 


Kitaev chain, for cr z =t, The susceptibility becomes: 

(A eff sin k) 2 p 


n(w) = -a 2 Yi 

k>0-,p,a=zL 


Eka U + ^PEka + ip ’ 


where are given by the Kitaev spectrum with: 


(15) 


nf a = n- yi 7sinfc) 2 + F| (7 , (16) 

A7 


A e ff = 


\J (7 sin fc) 2 + F| 


(17) 


All the results from the previous section apply to this case 
but with the ^-dependent parameters showed above. The 
system is in the topological nontrivial (trivial) regime for 
Vz > \J P 1 + A 2 (Vz < vV 2 + A 2 ). In Fig.|ijwe plot the 
imaginary part of the susceptibility as a function of the 
Zeeman splitting Vz for two different values of the cavity 
frequency uj. We see a similar behavior as in the case of 
the Kitaev chain: a peak emerges in U"(uj) at the topo¬ 
logical phase transition point, which becomes narrower as 
omega becomes smaller. However, an extra peak emerges 
at a larger Vz, and it is due to the resonance condition 
with the gaps around the k ~ kp (external gaps in the 
SO coupled nanowire spectrum). 


IV. MAJORANA FERMIONS DETECTION 

In this section, we consider a finite ID topological sys¬ 
tem coupled to the cavity (therefore with open boundary 
conditions, or OBCs), so that there are two Majorana 
fermions emerging in the topological region, each local¬ 
ized at one of the two ends of the chain. Taken together, 
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they give rise to a zero-energy fermionic state in the in¬ 
finite wire limit, which can be either empty or occupied, 
thus labeling the parity of a ID p-wave superconductor 
[18] . The Majorana wavefunctions decay exponentially 
in the wire on the scale of the superconducting correla¬ 
tion length £, and for a finite wire it can lead to a finite 
energy splitting cm oc exp (—L/£) of the initially zero en¬ 
ergy fermionic state pQ. In the following, we will show 
that both the presence of the Majorana fermions and the 
parity of the Majorana fermionic state can be inferred 
from the susceptibility II(cj). 

In the finite chain case we cannot obtain exact results 
for II(a;) anymore, therefore we proceed to calculate this 
quantity numerically (see Appendix [E]). We will treat the 
two models, the Kitaev chain and the SO coupled wire 
on equal footing, showing that they give similar results. 

For starters, the electronic Hamiltonian can be casted 
in the following form: 

H el = Ic+Mc, (18) 

with 

c = (KJ . . . {CiVsHcL} ... (cjvJ), (19) 

where s counts internal degrees of freedom, such as spin, 
band index, etc. For the Kitaev chain s = 1 (and thus 
we can disregard it), while for the SO coupled nanowire 
s =t,i- Here, M is a 2 Ns x 2 Ns matrix PQ, and we can 
write M = PWP 1 . with 

W 2p _ s , 2fe _ s = (-iy +1 6 p , k e k -, s = 0,1. (20) 

P is a unitary matrix ( PP T = P T P = 1) whose 
columns are the eigenvectors of M pQ. Also, e p , with 
p = 1,..., sN are the eigenenergies of the electronic 
Hamiltonian, including the Majoranas (if present). Thus, 
the electronic Hamiltonian can be re-written as 

H el = (1/2 )CWC, 
and C = P { C. where 

c = ({ci s } ... {cjvsHclj ... {ct a }), ( 21 ) 

with c£ (c p ) are the creation (annihilation) operators for 
the Bogoliubov quasiparticles in the finite wire, with p = 
1... N labeling the energy levels. Finally, we can write 

H e i = y^psC-ps ~ ^ j 5 ( 22 ) 

p,s ' ' 

and also define the spinorial wavefunction for the state of 
energy ±e ps at position j as tf ps (j) = (u J ps ,vl s ) T , where 
u ps( v ps) = H2j-i,p{P2j,p) are the electron (hole) compo¬ 
nents of the wavefunction at position j in the wire. 


The electron-cavity coupling Hamiltonian can be then 
written in the new basis as follows: 


H P i- ( 


= E 

p,p' 

x ( a t - 


nW rt c , / 

° ps,p's' c ps c p's' 


a ), 


' iG ps, p ’s’Clscl> s > + h -C- 


(23) 


(1 2 ) 

where Cp S ’ i s , are coefficients that depend on the trans¬ 
formation from the electronic basis C to the Bogoliubov 
basis C and read Hanoi: 


r (l, 2 ) 

^ps,p' s' 


N 

= a '^^ J ' l Pps(j) r z,y' l Pp' s' (j) • 
3 = 1 


(24) 


Here the pseudo-spin f = (r x ,r y , r z ) acts in the Nambu 
(or particle-hole) subspace. In general, all , s , 7 ^ 0 
for p p' and s 7 ^ s' , thus there are couplings be¬ 
tween all the levels (and bands) via the cavity field, 
and that includes transitions between the Majorana and 
the bulk (or gaped) modes. This in turn affects the 
correlation function in Eq. ([ 5 ]), which can be written 
as n(u;) = n BB (u) + U bm (uj) + n mmM, being the 
sum of the terms that contain only bulk states (bulk- 
bulk, or BB), cross terms between Majorana and the 
bulk (bulk-Majorana or BM), and Majorana contribu¬ 
tions only (Majorana-Majorana or MM), respectively. 
However, n mm(w) = 0 [32] due to the fact that the 
cavity cannot mix different parities, and in consequence 
the only contribution from the Majorana modes comes 
through the cross terms TLbm(w). We have found that 
for N 1 the U bb (uj) contribution is given by the one 
obtained from the PBCs in the first part of the paper, 
i.e., U BB oc N , while oc const , up to exponen¬ 

tially small terms in L/£. We note in passing that in a 
real wire, the smallness of the H B m compared toU BB is 
measured by Af/T, with X B being the Fermi wavelength 
and L the length of the wire. 

In the following, we analyze the cross-terms contribu¬ 
tion n bm{u). For cm e p with p 7 ^ M, we obtain: 


iismm = 

p,s^M 

x [|C (1) 12 


1 


1 


M,ps I ( n ^ n Ps) I P'M^ps 


e ps + id + if] e ps - id - irj 

(2) l2 (n M -l + rip S ) 


(25) 


where n ps and um are the occupations of the bulk and 
Majorana states, respectively. This is one of our main 
results. Inspecting the above expression, we see that it 
is strongly dependent on the Majorana state parity um • 
Assuming that e ps > 0 for p, s 7 ^ M and n ps = 0 for n 7 ^ 
M in the ground state, we obtain that H BM oc \C^ ps \ 2 

(IljjM \ C M,ps\ 2 ) for n M = 1 (n M = 0 ). To get more 
physical insight into the resulting susceptibility, we write 
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FIG. 4. Dependence of II'(a;,/x) on the chemical potential 
/x. The blue (dot-dashed), red (full), and the black (dashed) 
lines correspond to the susceptibility for PBCs, OBCs for 
parity um = 0, and OBCs for parity tim = 1, respec¬ 
tively. Lower Inset: the relative strength of the susceptibility 
All =1 2 |(n+ M — n BM )/(n^ M + n sM )| as a function of fi 
in logarithmic scale, for A = 0.1 (red-full), A = 0.2 (green- 
dashed), and A = 0.3 (blue-dot-dashed). The size of An' is 
exponentially reduced as a function of A. We used N — 50, 
uj = 0.2, A = 0.1, t = 1, and all energies are expressed in 
terms of t. 

(1 2 ) 

the coefficients C y Mps in the following way: 

^m] P s = + V M 2) Ups 

3 

— (u? M 5 r ,2 + v Mfir,l) V ps\ • ( 26 ) 

Let us analyze the implication of the above result. When 
e M = 0 , we also have u j M = and thus C$ p = C$ p , 
since electron and hole contributions are the same in the 
Majorana state. However, for a finite energy splitting 
cm 7 ^ 0 , and thus we have that u J M 7 ^ v J Ml which in 

turn results in 7 ^ (7^. All these suggest that the 
susceptibility n(cj), via n bm{u) should allow us to infer 
both the parity of the ground state and the zeros in the 
Majorana energy em, assuming their spatial overlap is 
large enough. 

In the main plot in Fig. [4j we show the real part 
n'( uS) for the Kitaev chain as a function of the chem¬ 
ical potential /x for the two parities um = 0,1 as well 
as the bulk value for PBCs. First of all, the values 
for n (uj) in case of periodic and OBCs are different be¬ 
cause of n bm(m), as this contribution has a different 
dependence on /i and A from the bulk states. Second of 
all, the open BCs wire susceptibility shows oscillations 
as a function of /i on top of the average value, of the 
form ±cos (kpL), with +(—) corresponding to tim = 1 
(um = 0 ): i- e - they are opposite in sign for the two par¬ 
ities. Here kp is the Fermi wavevector of the electronic 
system, and for the range of parameters considered is 
kp ~ 2/x pQ. This means that the cavity field can ac¬ 
cess the parity of the Majorana fermions non-invasively 
and without locally accessing the wire. Moreover, the 



FIG. 5. Dependence of n^cj, Vz) on the Zeeman splitting Vz- 
Above the topological transition (indicated in the figure), the 
red (full) and black (dashed) curves correspond to U'(u,Vz) 
for parity um = 0 and tim — 1, respectively. n'(cj, Vz) shows 
oscillations as a function of Vz, that are different in ampli¬ 
tude for the two parities tim = 0,1 (but having the same pe¬ 
riod), crossing at points where cm = 0 (see Inset). Below the 
topological transition the susceptibility reproduces well the 
one for PBCs. The parameters we took are [42] : L — 2/xm, 
t = 0.5 • 10 -2 eV, a = 0.4meV, /x = -l(T 2 eV, A = 0.25meV, 
uj — 0.02meV, and N = 80. 


oscillations disappear below the phase transition point 
\fi\ = 1, the susceptibility n^u;) acquires the same value 
as for the PBCs wire which signals that the Majorana 
fermions exist only above the topological phase transi¬ 
tion. In order to get a closer look at the oscillations of 
n(cj,/i), in the lower inset in Fig. [4] we show the real 
part of the relative difference between the two parities, 
An(w,/i) = 2 | (n^ M - n BM )/(n^ M + n BM ) | , for dif¬ 
ferent values of A. We see that the oscillations have 
the same periodicity as the Majorana energy splitting 
cm ~ exp (—L/£)| cos (kpL)\. Notice that the oscillations 
of the Majorana splitting with the chemical potential has 
been studied in detail [43| 44] together with the fact that 
the magnitude of the oscillations becomes exponentially 
suppressed in L/£ [45l 146] , 

In the main figure in Fig.[5j we plot the real part of the 
susceptibility for a ID topological wire as a function of 
the Zeeman splitting Vz for the two parities um = 0,1. 
The susceptibility n for that figure was computed us¬ 
ing realistic parameters that might be appropriate for 
an InSb wire such as in the experiments in Ref. |9j We 
find similar features as for the Kitaev toy model, namely 
oscillations as a function of the Zeeman splitting above 
the topological transition. These oscillations around the 
ground state have opposite sign and different amplitudes 
for each parity. Like for the Kitaev model, they have the 
same periodicity as the Majorana energy cm (see the in¬ 
set of Fig. [ 5 ]) and cross at points where cm = 0. Notice 
that if the parity is not conserved in the system (for ex¬ 
ample, due to the quasi-particle poisoning), n' will follow 
the ground state and exhibit therefore sharp cusps as a 
function A z at the crossing points where €m = 0 (see 
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FIG. 6. Dependence of II"(a;) on the cavity frequency uj for 
N — 50. The full-red (dashed-black) line corresponds to the 
susceptibility in the topological regime for g = —0.2, while 
the dot-dashed-blue (dotted-green) corresponds to the non- 
topological regime with open (periodic) BCs with /i — —1.8, 
so that the effective gap is the same A e ff = ||/i| — t| = 0.8 in 
both regions. Inset: a zoom in the region where the Majorana 
peak emerges. For all the plots we used t = A = 1, and all 
energies are expressed in terms of t. 


also Ref. EH for similar features in a topological Joseph- 
son junction). As for the Kitaev chain, we thus find that 
the cavity phase shift is thus able to detect the Majorana 
fermions and the parity of the ground state of a realistic 
topological wire. 

The imaginary part of II(a;) gives us also information 
on the presence of Majorana fermions. In Fig.[6j we show 
the dependence of II" (cj) on uj for the Kitaev chain, both 
in the topological and non-topological regimes, for t = A. 
We see that the Majorana fermions, through II bm(w), 
give rise to an extra peak in the susceptibility at half 
the effective superconducting gap A e ff = ||/x| — t| in the 
topological regime, while such a peak is absent for the 
same effective gap A e ff, but in the non-topological case. 
For completeness, we also show the result for PBCs, in 
which case there are no Majorana fermions. 

Finally, let us give some estimates for II(<j), and in 
particular for TLbm(w) and the resulting phase shift in the 
exiting photonic signal. We assume typical experimental 
values for the cavity frequency, uj c « 2 x 10 -5 eV, and 
with a quality factor Q « 10 5 , which results in photon 
escape rate k = 2 x 10“ 10 eV. For an estimate of the 
capacitive coupling a we refer, for example, to the case of 
carbon nanotubes, which have been under experimental 
scrutiny in the context of cavity QED |24[ [40]. There, it 
was found that a ~ 5.6 x 10 -7 eV, and we believe similar 
values should be relevant for semiconductor nanowires 
too. The phase shift of the radiation exiting the cavity 
satisfies 5(f) oc (a 2 /t k) so that we obtain 5(f) ~ 0.3 which 
is a sizeable value. 


V. CONCLUSIONS AND OUTLOOK 

We studied two paradigmatic examples of ID topolog¬ 
ical superconducting systems capacitively coupled to a 
microwave superconducting stripline cavity: the Kitaev 
chain and a ID nano wire with strong SO interaction in 
the presence of a magnetic field and in proximity of a su¬ 
perconductor. We analyzed the electronic charge suscep¬ 
tibility of these systems that is revealed in the photonic 
transport through the microwave cavity via its transmis¬ 
sion t(cj). We showed that this electronic susceptibil¬ 
ity can actually be used to detect the topological phase 
transition, the occurrence of Majorana fermions and the 
parity of the Majorana fermionic state in a non-invasive 
fashion. Such effects are due to the interplay between 
the bulk and Majorana states, either via virtual or real 
transitions taking place between the two, and which are 
mediated by the photonic field. As an outlook, it would 
be interesting to use the same cavity QED setup to ac¬ 
cess the physics associated with the fractional Josephson 
effect. 
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Appendix A: Derivation of the effective Kitaev 
Hamiltonian in the presence of the cavity field 

In this section, we provide theoretical arguments for 
the wire Hamiltonian utilized in Eq. ([6|, and the effec¬ 
tive electron-cavity Hamiltonian used in the Main Text 
(MT). In a continuum description, the natural way to 
account for the interaction between the electrons and 
the electromagnetic field is via the minimal coupling, 
i.e. p —> p — (e/c)A in the electronic Hamiltonian, with 
A being the electromagnetic field vector potential and 
p being the momentum of the electrons in the mate¬ 
rial. In a tight-binding picture instead, one accounts for 
the coupling between light and matter by performing the 
Peierls substitution to the hopping parameters ta+x be¬ 
tween neighboring sites i and i - 1-1, namely 

tii+1 tii+ dr , (Al) 

with A(r) being the electromagnetic field vector potential 
at position r, and the integration is performed between 
the sites i and i + 1. We will focus on the derivation of the 












effective Kitaev model in the tight-binding picture, as the 
microscopic, continuum model was described in great de¬ 
tail very recently in [40^. We thus refer the reader to that 
paper for a detailed calculation of the cavity effects, as 
well as the derivation of the capacitive coupling starting 
from the minimal coupling. 

Here we give some details on the derivation of Eq. © 
in the MT starting from a non-superconducting nanowire 
coupled to a bulk p -wave superconductor with such a cou¬ 
pling being assisted by the cavity field. For simplicity, we 
assume the bulk to be not s, but p-wave paired, thus the 
presence of spin-orbit coupling in the wire is not a nec¬ 
essary ingredient. However, the present calculations can 
be straightforwardly generalized to more realistic system, 
such as nanowires with SOI. The total Hamiltonian of the 
system reads: 


bulk and wire Hamiltonians. We get: 

H b = £fc,;>4.f> C fc,fr - E iAsinA: ( C ~M C M “ C k,b C -k,b) 

k k> 0 

(A6) 

H w = ^ ^ ^k,wc\' W Ck,w 5 (A7) 

k 

H t = Unt ^2 ( e ^4 ,w°k,b + e“^4 , b Ck, W ^j , (A8) 

k 

where = t p cos k — p p , with /i a the chemical potential 

in the p = w,b system. 

Next we perform the so called Lang-Firsov transfor¬ 
mation on the system Hamiltonian, which means H sys = 
exp (S)H sys exp(— S) with S chosen as follows: 


H sys — Hi) + H w + Ht + H c , (A2) 


S = a ] )'Y^c\ w c q ^. (A9) 

LU C q, 


where 


H p — p p ^2 c ],p c j,p + ^2 (tp c ],p c j+t,p 


^p c ],p c J+i,p ■ 




(A3) 


with p = b(bulk),w(wire ), and A w = 0 (no intrinsic 
superconductivity in the wire), and = A the p -wave 
pairing in the bulk superconductor. Here, Cj :P (ct ) and 
t p are the electronic annihilation (creation) operator at 
position j and the hopping parameter in system p = 6, re, 
respectively. The tunneling Hamiltonian in the presence 
of the cavity reads: 

Ht = Y, c] w c jib + h.c.) , (A4) 

j 


where (f)j = Ajdj , with Aj = i(aj /uj c )(at — a), dj, aq, cj c , 
and a (at), being the cavity vector potential, the cou¬ 
pling strength, the cavity frequency, and the cavity pho¬ 
ton annihilation (creation) operators, respectively. Note 
that we assumed that the cavity field points perpendic¬ 
ularly to the wire, and it has no component along it. If 
instead such components would exists, we should have 
modified the wire Hamiltonian too in order to account 
for the cavity induced phase factors. In the following, we 
will assume that ajdj = aj = <a, namely it is constant 
along the entire wire. Finally, the Hamiltonian of the 
cavity reads: 


H c = uo c °^ a > 


(A5) 


After some lengthy, but straightforward calculation we 
obtain the system Hamiltonian as follows: 


H sys H w T Hj) T a ^ ^ c^ w Cq^ w (a T u^) 

Q 

_ 

'V 

H c — w 

2 

T t{ n i + h - c -) +R c 

k 

" -V-' 

TV 2 Ht 

(A10) 

which implies we excluded the photonic field from 
the tunneling term at the expense of adding photon- 
dependent chemical potential shift in the wire (third 
term) as well as an interaction term (fourth term). Note 
that for ti n t = 0, the transformation does not affect the 
spectrum, as it can be simply undone. However, as will 
see in the following, in the presence of the tunneling term 
the photonic field in the form of the capacitive coupling 
can lead to real effects. 

In the following, we aim at finding an effective Hamil¬ 
tonian describing the wire only by integrating the bulk 
superconductor degrees of freedom up to second order in 
the tunneling ti nt . We choose to do so by employing the 
Schrieffer-Wolff transformation formalism, which means, 
as before, that we unitary rotate the system Hamiltonian 
as 

= e Ssw H sys e~ Ssw = H w + H b + H w _ c + H T + H c 
+ [5sw> H w + Hi + H w - C + Ht + H c \ + ... , 

(All) 



with uj c being the (fundamental) frequency of the cav¬ 
ity. Before deriving an effective wire Hamiltonian, it is 
instructive to switch to the Fourier space, for both the 


and choose 

[SswiH w + H h ] = -H T > (A12) 
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or Ssw = (£ w +£&) _ 1 #t, with C a being a superoperator 
whose action is defined as C a A = [H a , A ], \/A. This is 
equivalent to the following identity: 


Ssw = i lim 

77 —>-0 



(A13) 


This term excludes the tunneling Hamiltonian Ht in 
leading order (assuming there is no diagonal contribu¬ 
tion caused by such a term). Then, we neglect the con¬ 
tributions of the higher order terms on the wire spectrum 
by averaging over the bulk ground state | 0 &) in order to 
derive a purely (renormalized) wire Hamiltonian: 


Hf « (0 6 \H b +H W +H C + H w _ c + - [Ssw,H t ] 


+ [SW, H w _ c \ + ... |0ft), 


(A14) 


In order to find Ssw from Eq. ( |A13 ) explicitly, let us 
perform Bogoliubov transformation for the bulk p -wave 
superconductor defined as 


Ck,b = U*klk, b + Vk-y-kfi > ( Al5 ) 

C-k,b = ~ v klk,b + U k ^ ] _ k b , (A16) 

where k > 0 and 

Uk = \A/ 2 (1 + ik/Ekb) 

Vk = yi/2 (1 - Zk/Ektie-*' (A17) 


with (\>k the phase of the superconducting conden- 


sate (that we choose = 0 from now on) and Ekb = 
+ A 2 sin 2 k the spectrum. We can then express the 
bulk Hamiltonian in terms of the 7 k and 7 operators: 

Hb — Ekb {lk,b~Yk,b + • 

k> 0 V 7 

(Ais) 

Utilizing the fact that: 


Ck,w{t) = c k , w ( 0) exp (-i£k,wt ), 

(A19) 

and 


7 k,b(t) = 2fk,b (0) exp (~iE k ,bt ), 

(A 20 ) 


we can readily find the transformation matrix Ssw as 
follows (assuming also that £k,w £ 7 - 7 , since we are 
interested in the energies well inside the band gap of the 
bulk superconductor): 


Ssw = ^2 “ l^) ( C t,w C k,b - 4,ftC fc)U) ) 


^^kVk(Ck iW C—k,b C—k,bCk,w) 


(A 21 ) 


Utilizing this expression for Ssw, we can calculate 
the expec tation values for the different commutators in 
Eq. (A14). We obtain: 


H ind , w = hoftl [Ssw, Ht + 2 H c _ w ] | 0 b ) » - V §* 

1 k ^ kh 


X 


( \ u k \ 2 - M 2 ) c\ w c k}W - 2 UkV k cl k w c{ !W + h.c.j , 

(A22) 


which can be interpreted as follows: the first term renor¬ 
malizes the single particle spectrum in the wire, while 
the second term is responsible for the induced supercon¬ 
ductivity in the wire. The full wire Hamiltonian thus 
becomes: 


H w — ^ ' {£k,w H“ 8£k,w) c \ w C k, w 
k 

deff 
* k,w 

+ 2 E ( A i n d C k,w C -k,w + h ' C ) 

k 

a 2 

+ a J2 c k,w°k,w (« + « f ) + —Nl 


(A23) 


with 


Kk, w = (I n k \ 2 - \v k \ 2 ) = 

h kb Ekb 


t 2 

^ind = UkVk 

&kb 


1 2 
k/| nrt 


2 E 2 kb 


A sin k 


(A24) 

(A25) 


being the renormalization of the single-particle energies 
and the p -wave induc ed gap (oc sin k). Note that the 
last term in Eq. (A23) can be seen as a normalization of 
the single-particle spectrum in the mean-field, and thus 
finally we recover the same wire Hamiltonian defined in 
Eq. (| 6 | in the MT. 


Appendix B: Input-output theory for microwave 
cavities 


In this section, we present details on the input-output 
theory for the cavity in the presence of the coupling to a 
ID p -wave SC. We will show that the transmission of the 
cavity depends on the electronic susceptibility of the elec¬ 
tronic system. We concentrate on only single resonance 
of the cavity with frequency uo c [48| and, for simplicity, 
we assume for the moment a single-sided cavity. The 
total Hamiltonian describing the system reads 

H = H e i + H c + H e i_ c -\-Hfr + i^c-6 5 (Bl) 

v -V-' 

H S y S 
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C-out 



\AAAA/ 
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bout 


FIG. 7. Sketch of a one-sided cavity system probed in pho¬ 
tonic transport. The input fields bi n and Ci n are sent from 
the right and left mirrors, respectively, towards the cavity, 
and an output fields b ou t and c ou t are collected on the same 
sides. The cavity field, which is quantified by the bosonic 
operators a and a\ interacts with the electronic system via 
capacitive coupling, affecting the cavity. The coupling be¬ 
tween the cavity field and the external modes is quantified by 
the decay rate ^ 1 ( 2 ) = 27rp|/i( 2 )| 2 , with p the bath density 
of states and |/i( 2 )| 2 = f 1 ( 2 ) f 1 ( 2 ) being the average coupling 
between the cavity and the bath modes on the right (left). In 
the first part, for simplicity, we assume that ft 2 = 0, which 
correspond to a one-side cavity. 


where H e i is the electronic Hamiltonian only, and 


H c = uj c a)a , (B2) 

H ei - C = a(a + a t )n, (B3) 

H b = Y^ q b%, (B4) 

Q 

H c _ b = -ih Y, {f q a\ - f*b\a) , (B5) 

Q 


are the cavity, the electronic-cavity, the bath, and the 
cavity-bath Hamiltonians, respectively. Above, a (a^) is 
the annihilation (creation) operator for the cavity mode 
with energy cj c , b q (bj) are the annihilation (creation) op¬ 
erators for the bath modes with energy uj q , with q label¬ 
ing their quantum numbers, and the complex coefficients 
f q are coupling parameters between the cavity and the 
external bath. Moreover, a is the coupling strength be¬ 
tween the cavity field and the total number operator in 
the system n = c J c j> with c j (cj) being the anni¬ 

hilation (creation) operator for the electrons (fermionic 
degrees of freedom) at site j in the electronic system. 

The idea of the input-output theory is to find the out¬ 
put photons (or field) in terms of the input ones, as shown 
schematically in Fig. [ 7 ] Following Ref. [48], we obtain 
for the cavity equation of motion for a one-side cavity 


(« 2 = 0 in Fig. [7]) : 

d = — [Hgyg, a] —cl yj-ftbiri 5 (B6) 

for the input field, and 

d = — [Hgyg, cl] + ~cl \fftb ou t , (B7) 

for the output field, where k, = 27rp|/| 2 is the cavity decay 
rate, with p the bath density of states and / being the 
average coupling betw een the ca vity and the bath modes. 
Subtracting Eq. ( |B7| ) from Eq. ( |B6| ) we obtain that 

boutit) = b in (t) + V^a(t), (B8) 

a result which holds for any general cavity Hamiltonian. 

In the following, we will establish the relationship be¬ 
tween b out and bi n in the presence of the electronic sys¬ 
tem, as depicted in Fig.[7| For that, we first evaluate the 
commutator: 


[H sy3 ,a] = —LU c a - an , (B9) 

where n is the time-dependent electronic particle number 
(see below). In order to utilize this contribution to the 
equation of motion of the cavity field, we need to evalu¬ 
ate the time-dependent particle number n(£), which itself 
depends on the coupling to the cavity. At time £, we can 
write: 


nn(t) = e lHsys ( t to ^ne to \ 


and 


ni (t) = e lHel ^ to W lHel ( t to \ 


(BIO) 


(Bll) 


being the Heisenberg and interaction pictures, respec¬ 
tively, with t > to being the initial time and which can 
be chosen at will. We can represent n#(£) as follows: 

riH(t) = U\t,t 0 )n I (t)U(t,t 0 ), (B12) 

where 

U(t,t 0 ) = T c exp dt'H e i- c (t (B13) 


is the evolution operator with T c the time-ordering op¬ 
erator that puts operators with later times to the l eft of 
the ones with earlier times. We can then write Eq.(B12) 
in the following way 
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n H {t) = T c ex p ^ dt' H e i_ c (t')J m(t)T c ex p (-i^ dt'H e i_ c (t') 

« ( 1+iJt dt'H e i_ c (t')^n I (t) (l-iJt G?i / iy eJ _ c (i / )^ 

^n/(t)+z f dt' H e i- C (t r ),ni(t) = ni(t) + ia f dt' (a + (J)ni(t'), n/(t) 
Jtn •'to 


(B14) 


r 


up to leading order in the coupling constant a. Thus, ducing Eq. (B14) into Eq. (B9) we obtain: 
the time-evolution of the electronic particle number con¬ 
tains, besides the electronic component, a contribution 
that arises because of the coupling to the cavity. Intro- 


[. H sys , a] = — uj c a — an = — uj c a — ani(t ) 


ia 2 f dt' (a(t') + a\t')) n/(t / ), ni(t) . 

Jt 0 L J 

(B15) 


Let us assume that to — oo and write Eq. (B15) as: 


[H sys , a] = — uo c a — ani(t) — ia 2 f dt' (a{t)e luJc ^ ^ + a) {t)e luJc(<t t ^^/(t / ), n/(t) , (B16) 

J — oo L ' ' J 


where a(t) ~ ae~ luJct and a^(t) ~ a t e ^ c £ in zeroth order 
in a (because the expression is already multiplied by a 2 
we can utilize the bare time dependence in this expres¬ 
sion). In the following, we switch to the Fourier space in 


I 

order to solve the equation for a(t), and take into account 
all contributions that affect its time-dependence, namely 
the external modes too. We obtain: 


.a 


-i(jja(uj) = —iuj c a{ui) - a(uS) — ^/nbi n (uj) — i—ni(uj) 

a n 


f dte iut f dt^L^e-^'-V+a\t)e i “'( t '- t A ni (t , ),n I (t) 1. (B17) 

^ d — OO j — OO ' 


Before continuing with the derivation, let us describe 
each term in the above expression. The first term de¬ 
scribes the free cavity evolution, the second term the 
leaking into the continuum of modes (the external bath) 
at rate k/ 2, the third term is the input field supplied 
from the right side, the fourth term correspond to an¬ 
other “input” contribution to the cavity from the elec¬ 
tronic system (a noise term), while the last term leads 
to both a shift in the cavity frequency as well as to an 
extra decay channel (Q-factor modification). One can 
now average over the electronic system, thus neglecting 
any fluctuation (i.e. feed-back effects). Moreover, we can 
neglect the highly oscillating term a){t) oc e ZUJct , namely 
we perform the so called Rotating Wave Approximation 
(RWA). Under all these assumptions, the last term in 


r 


Eq. ( |B17| ) becomes: 

N 2 p°° 


a‘ 

h 


/ OO nt 

dte luJt / dt'a(t ) 
-oo J — oo 

a 2 r°° . f 

■t L***°®L 

([n/(0,n/(t)])o = i J 


— OO 
OO 


= ia(u)IL(— uj c ), 
where 


dt\-i)0{t’ - 
dte lwt a(t)H(—u> c ) 

(B18) 


u(t' - 1) = -iO{t’ -1) 


h 


{[n I {t’),n I {t)}) o, 


(B19) 
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is the retarded density-density electronic correlation 
function utilized in the main text, and (. •. )o means the 
expectation value of the unperturbed electronic system. 
Note that for deriving the above expression we assumed 
that the electronic system is in equilibrium and thus the 
time dependence is homogeneous. 

We are now in position to find the cavity field a(t) and 
the output field b out (t ) in terms of th e input field bi n {t). 
Introducing Eq. ( |B18| ) into Eq. ( |B17| ) we obtain 

Kj Ct 

—iuja(uj) = —iuj c a(uj) — —a(cj) — V^m(^) — 2 —(n/(cj)) 0 
+ uj c ) , (B20) 


so that 


V^bj n (uj) +i(a/H)(ni((jj))o 
— uj c ) + k/2 — zll(— uj c ) 


(B21) 


We thus have two contributions to the cavity field: the 
external input and the input from the electronic system. 
Our aim is to relate in fact the output and input fields, 
which can be done easily via the expression in Eq. (B8): 


bout (u) 

- cj c ) - k/2 - in(-u c )]b in (uj) - i(^a/h)(n I (uj))o 
—i(cj — uj c ) + k/2 — iU(—uj c ) 

(B22) 


In the limit of large number of photons in the input beam, 
we can neglect the contribution from the electronic sys¬ 
tem so that we obtain: 


bout (^) 


-i(u: - uj c ) - k/2 - iU(-u} c ) ( . 

•( \ 1 /q * tt / \Din\^) 

—i[uj — uj c ) + k/2 — Hl( — uj c ) 


(B23) 


In experiments, one actually encounters a two-sided 
cavity (see Fig. [ 7 ] for the nomenclature), in which case 
the expression for the cavity equation of motion reads: 


1 

a = — 


[Hsysi^l ^2 2 ) a \/ ^1 bin \/i 


so that for the output fields we get: 

(B24) 

bout — yjK\Q j bi n 

(B25) 

C-out — %/K 2 K + Ci n . 

(B26) 


By following the same reasoning as for the one-sided cav¬ 
ity, we obtain: 

/ \ _ V^l bi n (u) + y/K 2 Qn(^) + i(a/H) (nj(cj )) 0 

—-;--- 7 - 7 -——- 7 —, 

—i(uj — oj c ) + Ki /2 + K 2/2 — iTl(—oj c ) 

(B27) 


while if the two mirrors are the same Ki = k 2 = k , this 
becomes 

+ c mHl +i(a/h){n I (u>))o 
— i(u) — uj c ) + k — HI(—cj c ) 


a(cj) = — 


Assuming again that the input flux is much larger than 
the electronic contribution, we can write: 

c 0 utM = -Tb in (w) + (.. .)c in (u) (B29) 


with 

T - _ 1 _ - A e i<t> 

— uj c ) + k — iU(—uj c ) 


(B30) 


being the transmission of the cavity, which is a complex 
number, and which depends on the electronic suscepti¬ 
bility II(u;), as stated in Eq. 0 in the MT. 


Appendix C: The susceptibility for the Kitaev 
model in case of periodic boundary conditions 


In this section we give more details on the derivation 
of the susceptibility II (u ) for the case of a ring geom¬ 
etry for which we can apply PBCs. The Kitaev chain 
Hamiltonian in real space was defined in Eq. § as (for 
definitions of the parameters please see MT): 

Htot = H e i + H e i- C + H p h , (Cl) 

N 1 N—l 

Hel = -M E C i Ci - 0 E ( tC I C *+l + ACjCj +1 + h.c.), 
i=l i=1 

N 

H int =aJ24 c i( a + ai ), (C2) 

i= 1 


and H p h = uj c a)a. In this case, we can switch to the 
Fourier space, which implies we can write the fermionic 
operators as follows: 




(C3) 


where c& is the fermionic annihilation operator with mo¬ 
mentum k = 2irn/N. We can then rewrite the electronic 
Hamiltonian in momentum space 


H el = Y J H 1 4 dG {k), (C4) 

k> 0 


with 

H Bdo ( k ) = ik{c\c k - c- k cl k ) - iAsin h(c- k c k - c\c ] _ k ), 

(C5) 

is the Bogoliubov de Gennes Hamiltonian and = 
—t cos k — fjt. The interaction Hamiltonian between the 
electronic system and the cavity can as well be written 
in the k -space as 

Hei-c = E a (4 c fc _ c -kC ] _ k ){a + a f ). (C6) 

fc >0 

One can simply diagonalize the H e i in the k -space and 
write: 

H el = J2E k /rhk+4- k l-k), (C7) 

k> 0 
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with 


Ek = Eyj (—t cos k — n) 2 + (A sin k ) 2 . 
being the eigenenergies, and we used 
Cfc = WfcTfc +Vk1- k , 

c -k = -v* k 7 k + Ukl- k ■ 

Here, the functions Uk and vk are given by 


\ Vk \ 2 = -(i-E 

11 2 l E k 


(C8) 

(C9) 

(CIO) 

(Cll) 

(C12) 


We can now calculate the susceptibility U(uj) defined in 
the previous section, which quantifies the change in the 
photonic transmission r due to the interaction with the 
electronic system. In the Fourier space, we obtain: 


n (t) = -iO(t)a 2 ([ (4 c fc - c -fc c h 


( 


c r — ‘ 


k,q> 0 


J 


)(*)> 


,) (o)]). 


(C13) 


Usin g Eq . (C9) and Eq. ( C10| ) with Eq. (Cll) and 
Eq. ( C12| ) we find: 


x 9 (A sin k) 2 

n(w)~.c?Y - 2 >- 


k>0 


1 


1 


id — 2 El. 


■ IT] 


id + 2 Ek + if] 

( C14 ) 

with the small r] > 0 assuring the convergence of the 
time-integrals. For large N 1, we can transform the 
sum into integral, and also write n(cj) = n'(u;) +in // (cj), 
with: 


2Na 2 r 

nV) = / dk 

t Jo 


(A sin kY 


Eb 


to 2 - 4 El 


(C15) 


n» = 


N a 2 


f 


dk 


(A sin k) 


where V ... means the principal value of the function and 
we used the fact that: 


1 


x — a + ie 


= V- 


1 


x — a 


— i7rS(x — a). 


(C17) 


We can perform the integral over k for the imaginary part 
n'fcj) to obtain the expression in the MT: 


n» = 


a 2 tN . 
2 pad 


\ 


1 - -L 


(cj/ 2) 2 — t 2 — /i 2 


n 2 


2 ,, 2 


4t> 


(CIS) 


for |t + /i| < uj/2 < |t — /i|, and being zero otherwise. 
For the real part n '(id) we found no simple solution, and 
so we chose not to depict it. Note that the susceptibility 
n(cj) oc A, i. e. it scales linearly with the number of 
sites. 


Appendix D: The susceptibility of the nanowire in 
case of periodic boundary conditions 

The Hamiltonian for a nanowire in the presence of a 
Zeeman field H, Rashba spin-orbit coupling u and in¬ 
duced superconductivity A reads [45] 


H nw = (2 m~ ^) Tz+ up<JzTz + B(Tx + ^ Tx ' ( D1 ) 

Let us diagonalize the Hamiltonian Hnw = Ho + A r x 
in the absence of the induced pairing A and then treat 
the latter perturbatively 

Ho = - j Tz + nP^zTz + Ba x . (D2) 

In order to do so, let us perform the unitary transforma¬ 
tion 05] 


H 0 = UH 0 U\ 


(D3) 


where 


U = exp ( iaa y r z /2 ) = cos (a/2) + i(J y r z sin (a/2). 

(D4) 

If we choose a so that 

tan (a) = —, (D5) 

up 


Hq takes the from 


Ho = 


2 

^-/i + yju 2 p 2 + B 2 a z 

12 m 


(D6) 


Now let us reintroduce the pairing term A r x and apply 
the transformation U to it. 

Then Hnw reads 


[£(u; — 2 Ek) — S(cd + 2Ek)\ 
(C16) 


Hnw = 


—- p- 1- yju 2 p 2 + B 2 c 

12m 


t z + 


A up 


A B 


yj U 2 p 2 + B 2 


J y‘y • 


yj u 2 p 2 + B 2 

(D7) 


Neglecting the last term in Eq.(D7) the effective 
Hamiltonian reads 


H t 


eff 


^-/i + yj u 2 p 2 + B 2 cr z 

L2 m 


A up 


yj U 2 p 2 + B 2 


(D8) 


The bulk energy spectrum of H e q for a z = — 1 is 


e P,~ 1 = ±1 


- (I — yj U 2 p 2 + B 2 

2m 


A 2 2 2 

A up 

u 2 p 2 + B 2 
(D9) 
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and for o z = +1 is 


1. The projection onto the lower band cr z = — 1 


e P ,+1 =±\j [^- fi + p + -B 2 


A 2 2 2 

A u p 

uV + B 2 ' 
(DIO) 


The electronic susceptibility is defined as 

II(t) = -^(t)([r 2 (t),r 2 (0)]). (Dll) 


In order to diagonalize the Hamiltonian in Eq. (|D 8 j), 
let us perform a transformation 

H d = U 1 H eff Ut=e p T z , (D12) 

where 

U\ = exp ( i'yr y /2 ) = cos ( 7 / 2 ) + sin (q/ 2 ). (D13) 

Then the electronic susceptibility reads 

n(£) = -iO(t)( cos( 7 )r^( 0 ) - sin^r^t), 

cos( 7 )r 2 ( 0 ) - sin( 7 )r ;r ( 0 ) ). (D14) 

r x (t) = e iepTzt r x (0)e~ iepTzt 

= cos (2 e p t) r x ( 0) — sin (2 e p t) r y ( 0). (D15) 

Introducing Eq.( D15| ) into Eq.( |D14 ), U(t) reads 


n(t) = —2 6(t) sin 2 (7) sin(2e p £). (D16) 


Performing the Fourier transform, Eq. (D16) reads 


/* + oo 

n(cj) = -2 V sin 2 (7) / sin( 2 e p t) 

' „ Jo 


= - 2^2 sin 2 ( 7 ) 


= y]sin 2 (7) 


2e D 


4 e 2 + (7-M 2 


cj — 2 e p -\- ip uj + 2 e p + ip 


(D17) 


Sum can be transformed into integral as 

N 


5> £/*-£/**>• (D18) 


Then the imaginar y par t of the susceptibility n(cj) = 
n'(u;) -bin" (7;), Eq. (D17) reads 


The low-energy subspace at p = ±pp is formed by the 
bands for which a z = — 1 . pp = 2 rmx when p = 0 . 

The effective Hamiltonian for eso = mu 2 >> £>, p = 0 
linearized around ±pp reads 


H e ff = u(\p\ -p F )r z + sign(p)Ar x . 
The bulk energy spectrum reads 

e p = iy/A 2 + u 2 (|p| —Pf) 2 - 
In this limit 

sin 2 (7)(e) = 

and the density of states reads 

p(0 = 


u \/e 2 — A 2 sign (p) 


(D20) 

(D21) 

(D22) 

(D23) 


Introducing Eq.(D 22 ) and Eq.( |D23 ) into Eq.(D19) the 
imaginary part of the electronic susceptibility reads 


n» = 


N A 2 


V 'uj 2 — 4A 2 


(D24) 


urn 


2. The projection onto a z +1 

In this case the effective Hamiltonian linearized around 
p = 0 reads 


Heff = M + s) A + (D25) 

where an effective mass m* was introduced as 


1 _ 1 rP 
m* m + B 

The bulk energy spectrum reads 


(D26) 


e p *= [ B — p 


V 

2m * 


“bO 2 + c P 2 > 


A 2 u 2 p 2 

H 2 


where 


2,,2 


B — p A u 
- 7 -1~ 


n"M = J de P( e ) sin 2 (7)(e) [<*( e “ w / 2 ) “ 5 ( e + w / 2 )] 

(D19) 


m" B 2 
The density of state reads 

e 


p(0 = 


\J~c [e 2 — {B — p) 2 ] 


(D27) 


(D28) 


(D29) 
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And 


sm( 7 ) 2 (e) = 


A V [ e 2 _ {B _ m) 2] 

B 2 e 2 c 


(D30) 


Introducing Eq.(|D29| and Eq.(|D30j) into Eq.(|D19j) the 
imaginary part of the electronic susceptibility reads 


with ±e n being the eigenenergies of the BdG Hamiltonian 
and n— 1,... N. That pertains to the following diagonal 
Hamiltonian: 


N 

^ ^ Cm 
m=1 


C-m Cn 



(E7) 


n» = - 


NA 2 u 2 sfu 2 - 4 (B - /i ) 2 
4WE2C 3 / 2 


(D31) 


Appendix E: The correlation function for open 
boundary conditions 


where c m (cj^) are the annihilation (creation) operators 
for the Bogoliubov quasiparticles, which are defined as 
follows: 


c = P { c 


(E 8 ) 


In this part, we present details on the calculation of 
the susceptibility for a finite wire with open boundary 
conditions. In this case, the excitation spectrum changes 
compared to the previous case, as in the topological re¬ 
gion the zero energy Majorana fermions emerge. In the 
MT we found that the susceptibility reads: 


nn 



dte lu>t (\[ni(t), fij(0)]\) 


(El) 


and which can be written as 


and 


c (ci, , C 2 , C 2 , • • •, c/v, Gy) • (E9) 

It is instructive to introduce the wavefunctions (i) = 

(' u k ,i: Vk,i) , where u k ,i (v k ,i) = P 2 i-i,k (P 2 i,k), and which 
are describing the state k = 0,..., 2 N at position i = 
1,..., N in the lattice and accounts for the electron (u) 
and hole (v) components, respectively. That allows us to 
write: 


n(cj) = + n SM (o/), (E2) 

being the sum of a bulk susceptibility, that can be con¬ 

N 

Ci = ^ ^ ^2fc — l,i&k T 'U'2k,ic\ 

k=l 

, (E10) 

structed from only the bulk (or gaped) states, and cross 
terms that involve both bulk and Majorana states, re¬ 
spectively. We use the discrete lattice model to numer¬ 
ically diagonalize the Hamiltonian for an electronic sys¬ 
tem with N fermionic sites. This can be written in a 

N 

cj = ^ ^ ^2/c — l,iCk T ^2k,iCk 

(Ell) 

compact form as follows: 

k=1 



H,, = -c*Mc 


with 


c (ci, 64 , C 2 , C 2 , • • •, Cw, tyy-) 




J \T 


(E3) 


(E4) 


where M is a 2 N x 2 N matrix. Moreover, assuming all 
entries in the matrix are real, M is also symmetric, we 
can write it as follows: 


M = PWP\ 


(E5) 


where W is a diagonal matrix with eigenvalues on its 
diagonal and P is a unitary matrix (PP^ = P^ P = I ) 
whose columns are eigenvectors of M. The matrix W is 
ordered so that 


fe 1 0 

0 —€i 


W = 




0 

Vo 


e N 0 

0 —cat / 


(E 6 ) 


so that we can rewrite Eq. (El) in terms of Uk,i and Vk,i 
as follows: 


N 


n(cj) — ^ ^ (1 Tlk){^- ^rn)^2k—l,i^2rn—l,i 

i,j,k,m= 1 

^ (p ij 2k,j^2m,j V2k,j'U'2m,j) 

1 1 


K UJ + if] — 62 /c-l — t2m-l CJ + if] + 62 /c-l + C2m-1 

T (1 Clk)Tl m V2k—iP j 2k,j'^2m—l,j 

1 1 


+ ir] — e 2 k- 1 + e 2 m-i w + irj + e 2k -i - e 2m -1 

T ^/c(f ^2k — l,j^2m,j^) 

1 1 


^ + 62/c-l — ^2ra—1 CJ + IT] — 62 /c-l + 62 m -l 

T {p2k — l,j^2m— l,j ^2k— 1,l,j) 

1 1 


w ^ + e 2 fe-i + e 2m _i u + irj- e 2 k-i ~ e 2m -i 

(E12) 
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Next we extract from this expression only the Hbm{u) Let us now introduce the coefficients C^ s \ s = 1,2 de¬ 
component. This reads: fined in the MT: 

N N—l / , \ 


n bmh = e e 


i ~ 1 E + irj + e 2 fc-i -a; - irj + e 2 /c-i 

— (1 — rife — 71 m) V2k-l,i'U j 2M-l,i'U'2k,jV2M,j 

— V2h-l,i'U'2M-l,iV2k,j'U'2M,j + ,j V 2k,j 

— V2M-l,i'U'2k-l,i'V2M,j'U'2k,j 


— ( n M ~ Tlk) 


V2k-l,iU2M,iU2k,jV2M-l,j 


— V2k-l,i'U'2M,iV2k,j'U'2M-lJ + ^2M,i^2k-l,i^2M-l,j^2k,j 


’ V2M,i'U'2k-l,i'V2M-l,j'U'2k,j 


and, using that U2k-i,i = ^ 2 k,i, it can be simplified even 
further to give: 


N N -1 


n ,bmm = e e 


+ 


i 


hi=i fe=i 


— (jiM — 1 + n k ) 


u + ir) + e 2 fc-i -a; - irj + e 2 fc-i 


~ U2k,iV2M,iU2k,jV2M,j 


+ W2fe,^2M,i^2fe,jW2Mj “ ^2M,i^2k,i^2M,j ^2/c ,j 

+ U2M,iV2k,iV2M,jU2k,j 


+ { n M ~ Tlk) 


~ U2k,iU2M,iU2k,jU2M,j 


+ U2k,iU2M,iV2k,jV2M,j ~ ^2M^2k,i^2M,j^2k,j 

+ V2M,iV2k,i'U'2M,j'U'2k,j 


N 


(j( s ) _ y^(rt2M,A,l +^2M,i^,2)%,i 


i=l 


- {u2M,i^s,2 + ^2M,i^,l)^2fe,i , (E15) 


(E13) which we can utilize to rewrite IIbmO^) as follows: 


N-l 


n mH = E 


i 


i 


fe=i 


e 2 /c + CJ + if] e 2 k - U - iv 


(tim ~ rik) 


C (1) 


- (n M - 1 + n fe ) 


C® 


(E16) 


(E14) anc l which correspond to the expression Eq. (25) in the 
V ^"' / MT. 
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